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PREFACE 


As  part  of  their  continuing  program,  RAND’  s  Physical  and  Quantum 
Electronics  Group  selects  theoretical  topics  for  study  which  should 
be  of  ultimate  interest  to  the  Air  Force.  The  present  study  gives 
a  calculation  of  the  time  required  for  metals  to  vaporize  when 
subjected  to  a  laser  beam.  The  theoretical  study  presented  here 
should  serve  as  a  guide  to  future  experimental  work. 
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SUMMARY 


In  this  Memorandum  we  consider  tne  interaction  of  laser  light 
with  a  metal.  A  model  is  proposed  in  which  the  dominant  mechanism 
for  electron  de-excitation  in  transition  metals  is  by  the  emission 
of  optical  phonons ,  as  contrasted  to  the  simple  metals  where  the 
dominant  mechanism  is  by  the  emission  of  acoustic  phonons.  Cal¬ 
culations  based  on  this  model  give,  for  steel,  an  electron — optical 
phonon  relaxation  time  ~  10”^  to  lCT^  sec  compared  with  an  electron- - 
acoustic  phonon  relaxation  time  ~  lO-^  to  10*^  sec  for  a  simple 
metal  such  as  copper. 
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I*  imODUCTIC^ 

In  this  Memorandum  ve  investigate  a  mechanian  by  which  excited 
conduction  electrons  in  a  metal  may  interact  with  the  crystalline 
lattice.  This  interaction  causes  the  electron  to  transfer  its  energy 
to  the  lattice.  Specifically  we  will  be  interested  in  the  case  where 

all  of  the  conduction  electrons  per  unit  volume  are  excited  by  photons 

of  energy  of  the  order  of  a  few  ev.  For  n  (number  of  conduction  elec« 
22  23 

trons  per  cc)  ~  10  to  10  a  total  photon  energy  of  the  order  of 

3  k  * 

10  to  10  joules  per  cc  must  be  deposited  in  the  metal.  Energy  den¬ 
sities  of  this  order  of  magnitude  might  conceivably  be  supplied  by  a 
pulsed  laser  beam  (that  is,  an  intense  coherent  li^rt  pulse). 

What  is  ultimately  desired  would  be  to  describe  in  detail  the  in¬ 
teraction  of  a  pulsed  laser  beam  with  a  metal  upon  which  it  impinges. 
Broadly  speaking,  such  a  description,  if  complete,  would  include  a 
discussion  both  of  the  absorption  of  the  incident  radiation  by  the 
metal  and  the  response  of  the  metal  to  this  deposited  energy,  includ¬ 
ing  possible  damage  effects.  At  present  such  a  description  does  not 
seem  possible. 

In  the  present  work  we  separate  somewhat  artificially  the  proces¬ 
sion  of  events  into  absorption  (electron  excitation)  followed  by  elec¬ 
tron  de-excitation  and  lattice  excitation.  The  validity  of  this  ap¬ 
proach  depends  on  the  relative  magnitudes  of  the  laser  pulse  length 
(in  time)  and  the  electron  lattice  relaxation  time.  That  is,  as  photons 

For  example,  the  vaporization  energy  of  iron  is  ~  -  ev  per  atom 
(Ref.  1).  d 
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are  absorbed  and  excite  electrons  to  higher  states ,  these  electrons 
interact  with  the  ions  of  the  lattice  and  lose  their  energy* 

Grossly  this  is  manifested  in  a  localized  heating  of  the  metal;  melt¬ 
ing  or  vaporization  may  even  result*  Indeed,  if  such  change  of  state 
begins  to  occur  before  all  of  the  laser  energy  is  deposited,  the  lat¬ 
ter  portion  of  the  incident  radiation  encounters  a  dense  high- tempera¬ 
ture  vapor  (perhaps  a  plaana)  from  which  it  may  be  scattered  (re- 
radiating  energy),  as  well  as  absorbed*  Thus  this  portion  of  the 
laser  pulse  may  not  even  reach  the  solid  metal  to  excite  appreciable 
numbers  of  electrons*  Owing  to  the  extreme  complexity  of  the  complete 
process  we  make  the  somewhat  idealistic  assumption  that  the  laser  pulse 
is  deposited  (electrons  are  excited)  in  a  time  which  Is  short  compared 

to  the  electron- lattice  relaxation  time  (which  we  estimate  In  this 

* 

Memorandum  ). 

In  Section  II  we  give  the  assumptions  and  discuss  a  model  of  a 
metal  that  is  used  to  calculate  the  electron- lattice  relaxation  time. 

In  Section  III  we  compute  the  relaxation  time,  and  in  Section  IV  we 
discuss  the  results  and  give  sane  numerical  estimates  for  steel. 

* 

Some  early  crude  experiments  on  steel  indicate  that  our  assump¬ 
tion  may  be  reasonable* 
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II.  MODEL  FOR  ELECTRON  --  OPT I CAL- PHONON  INTERACTION 

O 

In  the  following  we  are  using  units  with  ti  »  c  =  1  and  e  1/137. 
We  will  denote  the  Rydberg  (meV2)  by  R. 

In  the  present  work  we  assume  that  all  of  the  conduction  electrons 

22  23 

per  cc  (~  10  to  10  )  in  a  metal  have  been  excited  by  absorbing  a 

photon  of  energy  of  the  order  (l/8)R  (corresponding  to  a  photon  wave 

o  * 

length  of  S'  7300  A) .  We  assume  that  each  electron  interacts  with  one 
photon  only.  (In  non -metallic  lattices — solids  with  band  separations 
large  compared  to  the  photon  energy — non-linear,  multi -photon  processes 
may  be  very  important  in.  determining  the  absorption  of  the  incident 
radiation;  this  will  not  be  investigated  here.)  Before  interacting 
with  the  photons  the  conduction  electrons  are  assrned  to  be  in  an 
equilibrium  Fermi  distribution.  That  is,  the  electrons  fill  a  sphere 
(in  momentum  space)  up  to  a  limiting  momentum  P^  (or  equivalently  to 
a  limiting  energy  E^). 

For  free  electrons  it  is  well  known  that  Ep  J^Pp  =  (2m  Ep)^2j 
is  given  by 

Ej,  =  (l/2m)(3*2n)2//3 

f  (1) 

e?  (l/8)R  for  n  —  1022  ^ 

At  this  point  we  wish  to  divide  various  metals  into  two  classes: 
transition  metals  (e.g.,  Fe,  Ni)  and  simple  metals  (e.g.,  Na,  Cu) . 
Transition  metals  are  distinguished  frcm  simple  metals  by  having  the 

— 

The  photons  will  not  appreciably  interact  with  the  ions  because 
their  velocities  are  much  smaller  than  the  electrons  (the  electromag¬ 
netic  coupling  is  proportional  to  the  current,  i.e.,  velocity). 
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d  shells  of  the  constituent  atoms  unfilled.  For  example,  the  copper 
atom  has  a  completely  filled  3d  shell  (10  electrons)  and  a  half -full 
4s  shell  (one  electron);  the  iron  atom  has  six  electrons  in  the  3d 
shell  and  a  completely  full  4s  shell  (two  electrons).  In  the  band 
theory  the  transition  metals  are  characterized  by  having  two  unfilled 
bands,  namely,  a  wide  low-density  s- electron  band  and  a  narrow  high- 
density  d-electron  band  (in  the  simple  metals  the  d-electron  band  is 
full).  In  iron,  for  example,  the  4s~band,  whose  energy  functions 
approximate  those  of  free  electrons,  overlaps  the  3d-band,  and  can 
take  electrons  from  it.  It  can  be  shown,  by  various  experiments,  that 
there  are  ^  o.l  electrons  per  atom  in  the  4s-band  in  iron  (See  Ref.  2, 
Chapter  II,  Table  2.10).  In  the  transition  metals,  the  Fermi  level  of 
each  band  is  at  the  same  height  in  the  equilibrium  state  (cf.  Fig.  1). 

Next,  consider  a  metal  during  the  time  when  the  photons  are  being 
deposited.  In  the  simple  metals  the  conduction  electrons  ( s-electrans) 
are  induced  to  make  intraband  transitions  by  the  laser  light  (oJlaser 
=  R/8  es  1.7  ev).  The  excited  s-electrans  interact  among  themselves 
through  their  mutual  coulomb  interaction.  The  s-electrans  can,  to  a 
fair  approximation,  be  treated  as  free  particles,  each  with  an  ef¬ 
fective  mass  as  mass  of  a  free  electron.  The  electron- electron  relaxa- 

po  po 

tion  time  in  a  free  electron  gas  (of  density  —  10  to  10  vcc)  is  of 
-13  -14 

the  order  10  to  10  sec.  Therefore,  for  laser  pulse  lengths 
-13 

»  10  sec  we  expect  that  the  s-eleetrons  will  be  adiabatically 

- 

The  d-bands  are  rather  narrow  in  energy  compared  to  s-bands  be¬ 
cause  the  corresponding  atomic  wave  functions  do  not  spread  far  out 
in  the  lattice  to  interact  strongly  with  neighbors. 
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s -  band  d-band 


Fig.l  — Filling  of  levels  in  transition  metals 


heated  (up  to  a  final  temperature  given  by  the  energy  of  the  laser 
pulse  divided  by  Boltzmann*  s  constant)  by  the  laser  photons*  That  is, 
the  excited  electrons  can  be  taken  to  be  in  a  Fermi  distribution  at 
a  temperature  given  by  the  energy  in  the  laser  pulse  divided  by  Boltz¬ 
mann  *  s  c  on  st ant . 

The  excited  electrons  may  now  interact  with  the  ionic  lattice, 
giving  up  their  energy  to  lattice  vibrations  (i.e.,  to  phonons).  In 
this  case  the  dominant  decay  mode  is  known  to  be  through  acoustic 
phonons  (see  Ref*  2,  Chapter  9).  Electron- -acoustic- phonon  relaxation 

times  are  of  the  order  10  ~  sec.  Hence  in  a  time  very  small  compared 

-8 

with  the  laser  pulse  length  (we  envision  laser  pulses  £  10  sec)  we 
may  assume  that  metallic  bonds  are  broken. 

Now  consider  the  transition  metals  (e.g.,  Fe).  Because  of  the  un¬ 
filled  d-band,  electrons  are  able,  through  their  counteraction  with 
laser  photons,  to  make  interband  transitions  from  the  s-baad  into  the 
d-band,  as  well  as  intraband  transitions  to  higher  states  in  the  s-band. 
Electrons  in  excited  s-band  states  may  interact  among  themselves  via 
their  mutual  coulomb  interaction.  However,  because  of  the  very  high 
density  of  states  in  the  d-band,  these  excited  s-band  electrons  will 


* 

The  energy  of  acoustic  phonons  is  <  =  Kmaxv>  where  K-max  is 

~  n/lattice  spacing  —  10*  cm,  and  v  is  the  speed  of  sound  in  the 

c 

metal  ~  5  x  lCr  cm/sec.  A  typical  value  for  co  is  ~  0.1  ev. 

^  max 

Generally  the  energies  of  the  acoustic  phonons  that  result  from 

excited  electrons  are  one  or  two  orders  of  magnitude  less  than  u)  (i.e., 

_  max 

they  have  energy  of  the  order  10"'3  to  10~^  ev).  Hence  electrons  must 
undergo  multiple  collisions  with  the  lattice  before  enough  energy  is 
supplied  to  the  lattice  vibrations  to  break  the  metallic  bands  (£  1  ev  per 

p  o 

atan).  Assuming  it  takes  10  to  10°  collisions,  and  treating  each  col¬ 
lision  as  an  independent  event,  we  see  that  the  time  required  for  the 
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scatter  into  the  d-band.  Once  an  electron  is  in  the  d-band  it  is 
trapped,  so  to  speak.  That  is,  in  the  d-band  the  mutual  coulcmb  inter¬ 
action  between  electrons  is  much  weaker  than  for  s-band  electrons. 

This  is  because  d-band  electrons  are  localized  around  the  icn-core, 
as  opposed  to  s-band  electrons  which  range  over  the  whole  unit  cell. 

So  once  an  electron  gets  into  a  d-band  it  does  not  interact  with 
other  electrons  appreciably.  Therefore,  we  assume  that  the  excited 
electron  configuration  in  the  transition  metals  consists  of  the  d-band 
filled  as  far  as  possible , and  essentially  no  electrons  in  the  s-band. 

In  this  case,  electron  de-excitation  via  acoustic  phonons  will  not 
take  place  to  any  great  extent.  This  is  because  if  excited  electrons 
decay  to  the  s-band  by  emitting  acoustic  phonons  (which  have  energies 
«  1  ev),  these  electrons  (l)  will  interact  among  themselves  via  their 
mutual  coulomb  interaction,  and  (2)  they  will  absorb  acoustic  phonons. 
Both  effects  thus  tend  to  scatter  the  electrons  back  into  the  high 
density  of  states  d-band. 

On  the  other  hand,  electron  de- excitation  via  optical  phonons 
will  result  in  much  larger  energy-  transfer  to  the  lattice  per  colli¬ 
sion  than  would  de-excitation  via  acoustic  phonons.  In  fact,  the 
energies  of  optical  phonons  are  an  the  order  of  an  ev.  That  is,  in 
a  single  collision  (with  the  lattice)  the  excited  electrons  are  able 


metallic  bands  to  be  broken  is  of  the  order  10~10  to  10”11  sec,  which 
is  still  small  compared  with  the  laser  pulse  length.  We  can  therefore 
conclude  that  via  acoustic  phonons  a  simple  metal  becomes  a  high- density, 
high-temperature  vapor  in  a  time  short  compared  to  the  laser  pulse 
length. 
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to  lose  all  of  the  energy  given  to  them  by  the  laser  photons.  In  this 
case  we  expect  that  after  a  time  of  the  order  of  the  electron --optical- 
phonon  relaxation  time,  the  laser  energy  will  appear  as  lattice  vibra¬ 
tional  energy. 

Summing  up,  we  are  assuming  that  in  the  transition  metals  the 
dominant  decay  mechanisms  for  excited  electrons  will  be  through  optical 
phonons  (as  contrasted  to  the  dominant  decay  mechanism  for  excited 
electrons  in  simple  metals,  which  is  through  acoustic  phonons).  Al¬ 
though  the  validity  of  this  assumption  can  be  investigated  in  more 
detail,  its  ultimate  verification  must  rest  with  experiment. 

Optical  phonon  transitions  are  the  important  mechanism  for 
electron  decay  in  dielectrics  where  energy  gaps  (of  the  order  an  ev) 
separate  the  ground  state  from  the  excited  state.  Although  the  formu¬ 
las  developed  below  apply  to  dielectric s,  as  well  as  transition  metals , 
we  shall  give  numerical  estimates  for  steel  only. 

For  the  optical  branch,  the  relation  between  energy  and  momentum 

is  go  =  constant  (i.e.,  approximately  independent  of  K),  We  take  the 

optical  branch  to  consist  of  a  band  of  levels  of  energy  of  the  order  of 

* 

an  ev.  The  group  velocity  of  an  optical  mode  is  very  nearly  zero. 

This  means  that  if  the  electrons  are  excited  in  a  localized  volume  the 
optical  phonons  which  result  via  electron  de-excitation  will  stay  con¬ 
fined  to  this  volume.  That  is,  there  will  be  essentially  no  transport 
of  lattice  energy  as  there  would  be  for  acoustic  inodes  (for  which  the 
group  velocity  equals  v) ,  Also,  w^e  consider  the  unit  cell  to  be  made 
up  of  two  atoms  (e.g.,  body-centered' cubic  structure). 

* 

The  center  of  mass  of  each  unit  cell  remains  fixed  for  optical 


inodes 
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III.  CALCULATION  CF  ELECTRON— PHONON  DECAY  TIME 


In  this  section  we  calculate  the  transition  rate  (probability 
per  unit  time)  for  the  excited  electrons  to  decay  by  emitting  a  phonon. 

The  hamiltonian  for  the  system,  including  electron- -phonon  inter¬ 
action  but  neglecting  electron— electron  and  ion— ion  interactions  is 

H  ■  v  +  Eelee  +  Hlnt  <2) 


where  is  the  hamiltonian  for  free  motion  of  the  lattice,  Helec 

is  the  hamiltonian  for  the  free  electrons  (i.e.,  Helec  describes  the 
motion  of  the  conduction  electrons  in  the  periodic  field  of  the  fixed 
lattice)  ^  and  in  the  rigid  ion  approximation,  is  given  by 


Hint 


-  2  Vfr-E-i  J-ST(r-L) 
L,b  X,D  L 


(3) 


In  Eq.  (3)  V(r  -  L)  is  the  potential  of  the  undistorted  lattice  at  r 
due  to  atans  in  the  L’th  unit  cel].,  and  b  is  the  displacement  from 
equilibrium  of  the  b'th  at  an  in  the  L'th  unit  cell.  In  reality,  when 
the  nucleus  moves  the  surrounding  electron  core  moves  less. 

This  fact  is  neglected  in  Eq.  (3) • 

Expanding  Eq.  (3)  in  powers  of  Q  we  find 


H. 


int 


=  -  S  Orb 

L,b 


V  V(r  -  L)  +  5 


2 

L,b 


(ql, 


b^or 


"We 


x  v  v„  v(f  -  E)  + 

a  p 


For  the  moment  we  drop  all  but  the  first  term  of  Eq.  (4).  Expanding 
^  in  terms  of  phonon  creation  and  destruction  operators  (denoted 
by  (A?)*  and  respectively)  and  substituting  into  Eq.  (4)  we  find 
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H 


int 


iH~Z 

L  jb  j<y 


-iK*L 

e 


+  complex 


conjugate]  e£  •  VV(r  -  L) 

XV  y  D 

where 


(6) 


is  the  phonon  wave  vector 

is  the  frequency  of  the  phonon  of  wave  vector  K 

is  the  phonon  polarization  direction 

is  a  unit  vector  in  the  direction  of  motion  of  the 
bfth  atom  for  a  phonon  of  wave  vector  K  and  polarization  a 

is  the  mass  of  the  atom  located  at  b 

The  first  and  second  terms  in  the  square  brackets  in  Eq.  (6)  correspond 

to  creation  and  absorption  of  a  phonon  of  momentum  K,  respectively. 

The  transition  probability  per  unit  time  for  electrons  of  momen¬ 
tum  p  decaying  to  a  state  of  momentum  p1  and  emitting  a  phonon  of  mo¬ 
mentum  K  =  (p  -  p1 )  is  given  by 


K 

'“K 

a 

%b 


r.  =  2*|m|2  d 


(7) 


where  D  is  the  density  of  final  states  per  unit  energy  interval.  For 
a  two-particle  final  state  D  is  given  by 


D  =  (2k)3  f 


a3?. 


3.! 


d3? 

r  *2 


(2*)J  -  (2k) 


(8) 


The  matrix  element  between  initial  and  final  states  of  H.  ,  is  denoted 

int 

by  M.  The  initial  state  is  an  electron  of  momentum  p;  and  the  final 
state  is  an  electron  of  momentum  p1  and  a  phonon  of  momentum  K.  Ex¬ 
tracting  from  (Eq.  6)  the  term  appropriate  to  creation  of  a  phonon 
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of  momentum  K  ( see  Fig*  2)  we  have 

\  1/2 


M  =  il 


L  j"b 


•‘“'L  lb'(JVrp*VV(?-£)'l'pd3r  (9) 


The  integral  in  Eq.  (9)  extends  over  unit  volume*  In  Eq.  (9)  the  Y  *s 

P 

are  the  well-known  Bloch  wave  functions  for  electrons  in  a  periodic 
lattice*  (Y  =  eip#r  u  (r),  where  u  (r  +  L)  «  u  (r)  and  !y  inte- 
grates  to  unity  over  unit  volume*)  In  Eq*  (9)  u^  corresponds  to  a 
d-band  wave  function  and  u^t  to  an  s-band  wave  function.  The  eigen- 
value  corresponding  to  ^  is  p  /2m,  where  m  is  an  effective  mass. 

That  is,  the  effect  of  the  rigid  lattice  on  the  conduction  electrons 
is  to  cause  a  shift  in  their  masses;  otherwise  the  electrons  behave 
as  if  they  were  free.  In  the  following  we  take  m  ~  free  electron 
mass  3*8  x  10^R) . 

Owing  to  the  periodicity  of  the  lattice,  the  integral  appearing 
in  Eq.  (9)  may  be  transformed  as  follows:  The  electronic  wave  func¬ 
tions  must  satisfy  Bloch* s  theorem,  i.e., 


y?  -  £)  =  e_ip’L  y?) 


Each  term  in  the  sm  over  L  in  Eq.  (9)  contains  an  integral  over  the 
whole  volume .  If  we  change  the  origin  of  the  variable  r  in  each  inte¬ 
gral  to  the  center  of  the  L’th  cell  and  then  apply  Bloch's  theorem  we 
have 


M  =  i  HZ 
L,b,a 


1/2 


^NMb 


+  K  -  P) 


X 


(J  »*,  W(f)  TP  d3r) 


(10) 
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Since  the  integral  must  be  the  same  in  any  cell,  the  sum  over  L  may 
now  be  done.  This  gives  us  a  factor  N  times  a  Kronecker  delta 
5(K  +  p1  +  5  p)  (G  is  an  arbitrary  reciprocal  lattice  vector,  i.e., 
L  *  G  «  integer  multiple  of  2jt)*  The  factor  N  may  be  absorbed  by 
assuming  that  the  electron  wave  functions  are  now  normalized  to  unity 
over  the  unit  cell  and  not  over  unit  volume  as  was  previously  done. 

We  have 


M«i 


5(k 


+  P’  +  G  -  p) 


n 


1/2 


a^NMb 


X 


V  V(r)  Yp 


(11) 


Henceforth  we  omit  the  Kronecker  delta,  keeping  in  mind  that  in  a 
transition  ve  must  have  p  =  p‘  +  K  +  G. 


Now  we  specialize  to  optical  phonons.  For  optical  phonons  ve  re¬ 
place  by  a  constant  Summing  over  b  ve  have  (e^  ^  g  »  e^) 

1^2 

M  =  1(1/2  »0  N)l/2  (M^1/2  -  M^1/2)  E  U  (p,  p* )  (12) 

or 

where 

Utf(p,  P1)  -  ^  •  J  V  v  V(r)  Yp  d3r  (13) 

and  and  are  the  masses  of  the  two  atoms  of  the  unit  cell* 

Fran  Eq.  (12)  we  see  that,  to  first  order  in  the  ionic  displace¬ 
ments,  if  «  Mg,  M  ^  0*  This  is  understandable,  since  for  optical 
modes  the  center  of  gravity  of  the  unit  cells  remains  fixed  and  thus 
the  contrary  motions  of  the  atoms  just  balance  out  in  their  scattering 


effect,.  By  taking  into  account  multi-phonon  processes  (i.e.,  the 
effects  of  higher-order  terms  in  the  electron  potential  energy  (see 
Eq.  4))  one  obtains  a  non-zero  contribution  to  M.  Multi-phonon  pro¬ 
cesses  are  discussed  in  the  Appendix.  For  impure  metals  like  steel 
or  alloys  the  right  hand  of  Eq.  (12)  is  nonzero. 

We  shall  now  study  the  problem  of  evaluating  the  integral  de¬ 
fined  by  Eq.  ( 13)  •  In  order  to  evaluate  we  make  the  assumption  that 
the  potential  field  surrounding  each  ion  is  rigidly  attached  to  it, 
and  moves  bodily  with  it.  That  is,  we  assume  that  V(r)  can  be 
written  as  a  sum  of  atomic  potentials 

V(r)  =  S  Vft(r  -  x±)  (14) 

centered  upon  the  ions  at  x^.  Then 

(15) 

The  integral  in  Eq.  (15)  should  be  carried  out  over  the  whole 
volume.  However,  it  should  cause  little  error  if  we  take  the  integral 
over  the  unit  cell.  The  potential  of  the  bare  ion  does  go  beyond  the 
unit  cell.  But  the  other  conduction  electrons  tend  to  screen  out  the 
ionic  potential  beyond  the  unit  cell. 


In  the  model  used  here  (i.e.,  rigid  ion  potential  with  spherical 
symmetry  round  each  atom)  M  s  0  to  first  order  (in  phonon  coordinates). 
In  order  to  obtain  a  nonzero  first-order  effect  it  is  necessary  to  use 
a  model  where  the  non- equivalence  of  the  fields  about  each  atom  in  the 
unit  cell  is  taken  into  account. 

Calculations  based  upon  a  less  restrictive  assumption  have  been 
made  by  several  authors,  in  particular  by  Bardeen  (Ref.  3)>  but  in  the 
present  state  of  the  theory  these  are  refinements  which  are  of  doubtful 
value  since  they  only  affect  the  numerical  values  of  parameters  which 
cannot  in  any  case  be  calculated  very  accurately. 
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We  are,  in  effect,  using  the  Wigner-Seitz  method^  (also  known 
as  the  cellular  method).  In  this  method  each  atcm  of  the  lattice  is 
surrounded  by  its  Wigner-Seitz  cell- -the  cell  bounded  by  the  perpen¬ 
dicular  bisectors  of  the  interatomic  vectors  (e.g.,  for  a  lattice  with 
the  body-centered  cubic  structure,  the  Wigner-Seitz  cell,  is  a  truncated 
octahedron).  The  electron  wave  functions  (for  a  given  band)  are  of  the 
form 


where 


¥  =  eip#r  u  (r) 

P  ° 


uQ(r  +  L)  =  uQ(r) 


and  satisfies  the  boundary  condition 


(16)* 


duo/dn  =  0, 

S/Sn  denoting  differentiation  normal  to  the  bounding  planes.  Since, 
however,  the  Wigner-Seitz  cells  approximate  closely  to  spheres,  it  is 
a  good  approximation  to  apply  this  boundary  condition  over  the  surface 
of  a  sphere  of  equal  voltcne  (if  the  radius  of  such  a  sphere  is  r  ,  we 
have  4:tr':/3  =  atomic  voltxne) .  The  boundary  condition  thus1  becomes 


"3F 


-  0 

r*r 

o 


(17) 


Integrating  the  right-hand  side  of  Eq.  (15)  by  parts  we  obtain 

Ua(p,  p')  -  *  J*  va(?)  V  (Tp,  y  <*3r  (18) 


The  electron  wave  functions  for  different  bands  have  the  quali¬ 
tative  form  of  Eqs.  ( l6) .  However,  they  differ  quantitatively  by  hav¬ 
ing  a  different  function  u  (r)  for  each  different  band. 
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where,  now,  the  integral  is  taken  over  a  single  cell.  The  integrated 
terra  vanishes  because  the  integrand  is  periodic.  Now 

V  Yp  =  ip  Yp  +  eip*r  VuQ 

=  iptp  (19) 

where  it  is  assumed  that,  over  the  greater  part  of  the  volume  of  the 
cell,VuG  negligible.  Hence  from  Eq.  (l8)  and  (19)  we  obtain 

Ua(p,  p’)  a  ieK  •  (p  -  p*)  J  Va(r)  ¥*,  Yp  d3r  (20a) 

cc 

=  ieK  •  (R  +  G)  J  va(r)  **,  Yp  d3r  (20b) 

where  we  have  used  conservation  of  qua  si -momentum  in  obtaining  Eq.  (20b) 
from  Eq.  (20a)* 

Now,  Unp-Klapp  processes  (i.e*,  processes  for  which  G  ^  0)  are 
responsible  for  scattering  of  electrons  through  angles  greater  than 
2  sin  ^(2  2^)  as  790 #  (5)  jjere,  electrons  will  be  scattered  through 

angles  less  than  1  rad.  (~  57°)*  Henceforth  we  take  G  =  0.  Conser¬ 
vation  of  momentum  now  reads 

p  =  K  +  p1  (21) 

With  G  =  0  the  initial  factor  in  Eq.  (20b)  becomes  *  K.  This  indi- 

* 

cates  that  only  longitudinal  phonons  will  interact  with  electrons. 

In  order  to  estimate  the  integral  in  Eq.  (20b)  we  argue  as  follows; 
If,  for  the  moment,  ve  take  p  =  p1  (i.e.,  consider  an  intraband  transition 
with  K  =  0)  then  we  require  the  mean  value  of  the  potential  energy  over 

*  . 

For  processes  with  G  f  0  transverse  phonons  may  interact  with 
electrons.  For  then  there  is  no  guarantee  that  K  +  G  is  parallel  to  K. 
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the  unit  cell.  The  magnitude  of  the  mean  value  of  the  potential  energy 

is  of  the  order  of  the  average  electron  kinetic  energy,  which  is  of 

the  order  of  the  Fermi  energy.  On  the  other  hand,  consider  the  value  of 

the  integral  for  large  K.  Due  to  shielding  (by  other  conduction 

electrons)  V  (r)  a;  0  for  values  of  r  greater  than  the  order  of  the 
a 

interatomic  distance.  Hence  the  integral  vanishes  (because  the  inte- 

ilC .  r 

grand  contains  the  function  e  which  oscillates  rapidly  for  large 

K). 

Here,  Kr  £  1  which  is  intermediate  between  very  small  and  very 

large  K.  Also,  we  require  the  value  of  the  integral  for  describing 

a  d~band,  and  for  Y  t  describing  a  s-band  (i.e.,  we  want  the  value  of 
]? 

the  integral  for  an  interband  transition)  *  In  this  case,  the  initial 
and  final  wave  functions  contain  different  angular  dependencies,  which, 
on  integration  will  tend  to  reduce  the  value  of  the  integral.  Hence, 
we  expect  that  the  value  of  the  integral  is  less  than  E^  (perhaps  a 
fraction  <  1  of  E^) .  To  determine  exactly  how  much  less  we  would 
have  to  determine  a  realistic  potential  V(r)  that  describes  correctly 
the  properties  of  the  metal,  substitute  this  potential  into  a  Schro- 
dinger  equation  that  includes  electron- electron  interactions  ( so  as 
to  include  screening),  determine  the  electron  wave  functions  and  then 
substitute  the  potential  and  wave  functions  into  the  integral.  This, 
however.  Is  beyond  the  scope  of  the  present  work.  Furthermore,  since 
Kr  <  1,  we  see  no  reason  to  assume  that  the  integral  is  a  very  sensitive 
function  of  K.  Therefore,  we  take  the  value  of  this  integral  to  be 
independent  of  K,  and  of  magnitude  approximately  E^  =  R/8. 

Fran  what  has  just  been  said,  we  have 
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P1 )  =  iKI  for  longitudinal  modes 
=  0  for  transverse  modes 

where  I  a  r/8. 

From  Eqs.  (22),  (12),  (8),  and  (7)  we  have 


(22) 


T1  =  2jt  |M|2  D 

.  -  M^/2Hl2/2  »0»)  J  K2  5^ed(p)  -  ^ 


where  is  the  phonon  energy* 

In  arriving  at  Eq.  (23)  we  have  used  conservation  of  momentum 

(Eq.  (21))  and  integrated  over  the  final  electron  momentum  (using  the 

delta  function  of  momentum  that  appears  in  Eq.  (8))*  Here,  all  of  the 

2 

excited  d-electrons  (of  energy  c^(p)  =  +  V^,  where  is  the 

average  potential  energy  of  a  d-electron)  are  taken  to  he  in  a  narrow 

region  beginning  at  =  R/8.  The  final  energy  is  an  s-band  electron 

(of  energy  (p  -  K)  /2m),  and  an  optical  phonon  of  energy  u)Q  (which  is 

independent  of  K) .  Choosing  the  polar  axis  along  p  and  then  integrate 

* 

ing  over  the  azimuthal  angle  ve  find 


2  1 

ri  =  ^2)  (I2A«  u>0N)  J  dx  |  K4dK  S(p.ic/m  -  (K2/2m)  -  u)  )  (24) 


where  x  =  cosine  of  the  angle  between  p  and  K,  and  u)  =  u)Q-f  75^  -  e^(p). 


Integrating  over  K  we  find 


Fi  -  m^1/2)2(i2A«  »0B)  P3  J1  s  % 


^  /~2  2 
+  Vx  -  a 


a  (x*  -  a  ) 

(  /2  2~'4 
+  U  -  V  x  -  a 


175 


Hereafter  we  drop  the  subscript  s  on  m. 


(23) 


(25) 
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vhere 


2 

a 


=  2m  co/p 


2 


(26) 


The  lover  limit  on  the  integral  in  Eq.  (25)  comes  from  the  requirement 
that  K  be  real. 

The  integration  over  x  is  elementary.  The  final  result  is 

,  o  1/2  p 

ri  =  C  p3(l  -  a2)  (1  -  a  /2)  (27) 

where 

C  =  (M^1/2  -  M2-1/2)  (mI2/jtui0N)  (28) 


Now, 


10  "  ®o  +  S  '  ed(p) 


+ 


=  <n 


-  V, 


Substituting  Eqs.  (28)  and  (29)  into  Eq.  (27)  we  find 


(29) 


F,  sa  CP  3  (30) 

1  .b 

In  obtaining  Eq.  ( 30)  we  have  taken  u)Q  —  Vd  ~  Ep,  and  p  -  P^.  In  the 
next  section  we  discuss  Eq.  (30),  giving  sane  numerical  estimates  for 


steel. 
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IV.  DISCUSSION  AND  CONCLUSIONS 


Let  us  apply  Eq.  (30)  to  stainless  steel  (=  90  per  cent  Fe,  10  per 

cent  Cr).  In  this  case  we  assume  that  20  per  cent  of  the  unit  cells 

/£\* 

contain  a  Cr  atom.  Therefore,  the  interaction  amplitude  (Eq.  (ll)) 

should  contain  a  factor  2  X  IcT^"  and  the  transition  probability  a 
factor  h  x  1(T  .  Upon  substituting  the  Fe  and  Cr  mass  into  Eq.  (30) 
we  have  (we  take  =  l/8  R,  2N  ~  10^  cells/cc,  m  =  3*8  x  lO^R, 

I  «  Ep  =  R/8) 

^  =  l/Tl  5:  2  X  10'10R  (31) 


Now, 


R-2X  10l6 


-1 

sec 


Hence, 


2.5  x  10”^  sec  (stainless  steel) 


(32) 


For  a  steel  that  is  predominantly  90  per  cent  Fe  and  10  per  cent  Ni 
we  find 

/T  ^  -ft, 

T,  sb  1.5  x  10" D  sec  (Fe-Ni  steel) 


Physically,  in  time  an  excited  electron  has  a  probability  of 

=  37  per  cent  (i.e.,  e-1)  of  decaying  by  emitting  a  phonon  (i.e.,  the 

-  rit 

square  of  the  electron  wave  function  contains  a  factor  e  ).  In  a 
few  times  -r.  all  of  the  electrons  should  have  decayed  and  hence  the 


Cells  (80  per  cent)  that  contain  two  Fe  atoms  contribute  zero 
to  the  amplitude  (see  Eq,  (28)). 

—8 

For  a  50-50  Fe-Cr  alloy  t,  ~  1  X  10  sec ;  for  a  50-50  Fe-Mo 
alloy  T-^  ~  lo- 10  sec.  ~ 
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lattice  gains  energy  of  the  order  of  the  laser  energy  (here  taken 
3  4  .  * 

~  10  to  10  joules  per  cc).  Hence,  at  the  end  of  this  time  the 
lattice  energy  will  be  greater  than  the  vaporization  energy,  which 
implies  that  the  metallic  lattice  has  became  a  fluid.  This  fluid 
being  at  a  very  high  density  and  temperature  (  —  10^  to  10^°K)  will 
expand,  causing  large  pressure  gradients  (i.e.,  shocks)  etc. 

We  should  like  to  point  out  that  the  time  obtained  here  is  much 
longer  than  the  time  it  would  take  if  the  electrons  interact  with 
photons  whose  energy  »  R  (e.g.,  X-rays).  In  this  case  the  wave¬ 
length  of  the  excited  electron  is  much  smaller  than  the  interatomic 

0  1 

spacing  and  the  velocity  of  the  electrons  ~  10  cm  sec  .  The  elec¬ 
tron  makes  collisions  with  atoms,  giving  up  energy  to  the  atcms. 

Since  for  iron  the  vaporization  energy  (per  atom)  is  ~  l/2  ev,  it  will 
take  few  (say,  10)  collisions  (per  electron)  before  the  metallic  bonds 
are  broken.  The  time  in  which  this  happens  is  of  the  order  of  the  mean 
free  path  (~  10~^  cm)**  divided  by  the  velocity  of  the  electron 


'X' 

In  the  Appendix  we  show  that  the  transition  probability  for  an 
electron  decaying  by  emitting  two  phonons  is,  for  stainless  steel,  an 
order  of  magnitude  less  than  for  single  phonon  emission.  However,  for 
Fe-Ni  steel,  the  first  and  second  order  transition  probabilities  are  of 
the  same  order  of  magnitude.  This  result  is  quite  accidental  since 
in  the  first-order  transition  probability  for  steel  there  occurs  the 
difference  of  two  nearly  equal  numbers  and  a  small  factor  (4  x  10“2) 
which  takes  into  account  that  20  per  cent  of  the  unit  cells  give  a 
non- zero  scattering  effect.  The  first-order  transition  probability 
for  most  alloys  is  typically  much  larger  (decay  time  much  shorter) 
than  for  steel  (e.g.,  for  a  50-50  Fe-Mo  alloy  Tp  ~  10-10  sec)  and 
hence  the  first-order  scattering  dominates  (second-order  decay  times 
are  typically  ~  10"^  sec). 

**The  mean  free  path  =  (nc)"^,  where  the  absorption  cross-section 
a,  for  Fe,  can  be  taken  (at  X-ray  energies)  to  be  the  geometric  cross- 
section  of  the  Fe  atom  (a  —  10~l6  cm^).  Hence,  for  n  ~  10^2,  (na)"^*  ~ 
1CT6  cm. 
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g 

(  —  10  cm/sec).  That  is,  this  time  (for  ten  collisions)  is  of  the 
-15 

order  of  10  sec. 

Also,  we  note  the  manner  in  which  a  metal  "breaks  up  is  similiar  to 
the  X-ray  case  when  the  laser  photon  density  (per  electron)  is  much 
greater  than  one.  For,  when  the  photon  density  (per  electron)  is  much 
larger  than  1,  the  electron  picks  up  (on  the  average)  energy  much  larger 
than  the  Fermi-energy  (i.e.,  in  this  case  multiphoton  absorption  takes 
place)  or  equivalently,  momentum  much  larger  than  the  Fermi -momentum. 

That  is,  in  this  case  the  wavelengths  of  the  excited  electrons  are  much 
less  than  the  interatomic  spacing.  These  conditions  are  the  same  as  the 
X-ray  case. 

We  restate  that  optical-phonon  decay  is  the  principal  mechanism 
by  which  substances  with  band  gaps  (£  wXaser^  ^e“exc^e  dielectrics 

and  insulators).  Where  the  gap  spacing  is  much  larger  than  the  laser 
photon  energy,  the  material  will  be  transparent  to  the  laser  beam  un¬ 
less  the  average  energy  given  to  the  valence  electrons  is  £  gap  spac¬ 
ing  (i.e.,  we  would  expect  the  material  to  be  transparent  unless  the 
laser  photon  density  (per  electron)  is  large  compared  to  1,  so  that 
the  electron  may  undergo  multiple  photon  absorption).  For  further 
information  on  multi-photon  absorption  see  Ref.  7* 

We  emphasize  that  here  ve  have  calculated  the  redistribution  of 
excited  electrons  with  accompanying  phonon  excitation  under  highly 
non- equilibrium  conditions.  Only  when  the  photon  energy  is  deposited 

in  a  time  which  is  longer  than  the  electron-phonon  decay  time  (~  10~^  to 

-6 

10‘  sec  for  steel)  is  it  meaningful  to  think  of  the  metal  as  being 

vaporized  instantaneously.  Hence,  for  steel  and  with  laser  pulse  lengths 

-7 

<  10  sec  it  is  not  possible  to  apply  immediately  an  equilibrium  equa¬ 
tion  of  state  to  determine  pressure  gradients,  etc.  It  is  necessary  to 
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wait  a  time  of  the  order  of  the  electron-phonon  decay  time  before  such 

calculations  become  meaningful. 

Finally,  we  re-enrphasize  that  our  conclusions  about  the  decay 

mechanism  of  transition  metals  depend  completely  on  the  validity  of 

the  assumption  of  the  "d-band  trap"  (i.e.,  the  assumption  that  the 

conduction  electrons  become  inhibited  from  decaying  via  acoustic- 

phonons  because  of  their  affinity  to  stay  in  the  high-density  of 

states  d-band).  If  this  assumption  is  valid  then  we  conclude  that 

in  the  transition  metals  the  electron  de-excitation  will  take  place 

rj  6 

(via  optical-phonons)  in  a  time  ~  10  *  to  10  sec  as  contrasted  to 
the  simple  metals  where  electron  de-excitation  takes  place  (via 

acoustic  phonon)  in  a  time  ~  10  ^  to  10  ^  sec.  With  laser  pulse 

-8 

lengths  ~  10  sec  and  fast  photo- detectors  (with  resolving  time 
~  10~9  sec)  it  is  conceivable  that  experimental  investigation  of  the 
early  time  hi  story  of  the  metallic  decomposition  in  the  transition 
metals  may  be  done. 
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Appendix 
TWO-PHONON  DECAY 


In  this  appendix  ve  estimate  the  decay  rate  into  two  phonons* 
Two-phonon  final  states  come  about  in  two  different  ways  ( see  Figs.  3 
and  4). 

First  consider  the  process  depicted  in  Fig.  2.  The  amplitude 
for  this  process  arises  out  of  the  second-order  expansion  of  the  elec¬ 
tron  potential  energy.  Now, 


r2  '  2” 


m 


where,  frcm  Eq..  (4) , 

r(2)  _  1 


Hint  =  2  (^b^  V<*  V|3  V^r  "  ^ 


(35) 


and 


d%  d^P  d^P  ^  ** 

D  =  (2jt)3  [  - K  f  - 4r  f - \  &3  P -  -  P  .  ^  6  C  E,  -  E  \  (36) 

J  (io3  J  J  (io3  v  in  out^  v 111  out; 

Expanding  the  Q!  s  in  terms  of  phonon  creation  and  annihilation  opera¬ 
tors  ,  and  retaining  those  terms  that  correspond  to  the  creation  of  two 
phonons,  we  find 


H  p  '  ”  2  ' 


1/2 


/ _ L 


1/2 


ty  -O' 

^h  V,v 


J  tpl*7  W(r)  Yp  d3r 


X  8(p  -  p'  -  K  -  K  -  G) 


(37) 


We  denote  (HlQt)  by 

**Here  we  are  dealing  with  a  three-particle  final  state. 
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Fig. 3 —  Doubl  e  -  phonon  decay 


Fig. 4  — Two  successive  virtual  single-phonon  decays 
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In  obtaining  Eq..  (37)  we  have  performed  the  summation  over  L  (see  the 
derivation  of  Eq..  (10)  for  details).  Assuming  a  representation  for 
V(r)  such  as  in  Eq.  ( 14)  we  have  to  evaluate  integrals  of  the  form: 

J  VVVa(r)  Yp  d3r  (38) 

Where  the  integration  is  performed  over  only  the  unit  cell.  Inte¬ 
grating  Expression  (38)  twice  by  parts  we  obtain 

•** 

(K  -  K’)  (K  -  K’)  J  V(r)  T*,  Yp  d3r  (39) 

The  integrated  terms  are  zero  because  the  integrands  are  periodic 
over  the  lattice.  The  integral  appearing  in  Expression  (39)  is  identi¬ 
cal  to  the  integral  I  discussed  in  Section  III  (see  discussion  preced¬ 
ing  Eq.  (22)). 

For  optical  modes  Eq.  (37)  yields  (on  performing  the  summation 
over  b) 

rq 2  =  ~  u,0)(1At1  +  S  ^  *  (K  -  i')  *  (K  -  K* )  (4o) 

or 

Sane  qualitative  features  of  the  second-order  amplitude  may  be  seen 
from  Eq.  (39)-  First,  if  we  see  that  rr^  ^  0.  We  have  seen 

that  the  first-order  term  vanishes  if  (see  the  discussion 

following  Eq.  ( 13) )  •  Here,  the  ionic  displacements  appear  quadrat ically 
and  thus  the  scattering  amplitude  fran  each  ion  in  the  unit  cell  ha6 

-X- 

The  electron  wave  functions  are  assumed  to  be  normalized  to  unity 
in  the  unit  cell. 

As  before,  we  neglect  processes  where  5^0. 

***We  have  arbitrarily  set  tu^,  =  io^  =  u>0/2  where  u)q  ~  (l/8)R. 
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the  same  relative  sign.  In  pure  metallic  crystals  (e.g.,  pure  Fe) 
the  first  nan- vanishing  contribution  to  the  decay  rate  comes  from  the 
second-order  term  in  the  expansion  of  the  electron  potential  energy. 

In  order  to  compare  this  second-order  rate  with  the  first-order  rate 
for  steel  ve  put  **  Mg  =  mass  of  the  Fe  atom  (since  Fe  is  the  major 
constituent  of  steel).  Second,  ve  see  that,  in  the  second  order, 
transverse  phonons  will  couple.  Here,  however,  since  ve  are  interested 
only  in  an  order  of  magnitude  estimate,  ve  shall  neglect  the  contribu¬ 
tion  to  the  decay  rate  arising  from  transverse  phonons.  With  these 
simplifications  Eq_.  (4o)  becomes 

nr)2  =  -  (i/NM  u>0)  (K‘(K  +  K')/k)(K'  *(K  +  K*  )/K* )  (^l) 

Upon  substituting  Eqs.  (36)  and  (4l)  into  Eq.  (3*0  we  find 

Tg  =  (I/NM  tuo)2  J  d\  d^'/C2*)5  {s(ed(p)  -  [(p  -  K  -  K')2/2mj-  u)Q) 

X  (K  •  (K  +  K*  )/K)2  (K*  -  (K  +  K')A’)2}  (^2) 

vhere  the  upper  limit  on  these  integrals  is  —  K^Ry  (Kmax  ~  it/lattice 

spacing).  Changing  from  the  variables  K,  K' ,  to  K,  §(  =  K  +  K’ )  and 

-  -  2 

using  conservation  of  mcmentvm  to  express  K'  *Q  and  K'  in  terns  of 
Q  and  K  ve  find 

r2  «=  [l/(2jt)5]  (I/NM  a»Q)2  J  <13Q  6(p.5/m  “  Q2/2®  -  <») 

X  J  d^Q2  -  Q  •  ic)2(Q  •  K)2/K2(Q2  +  K2  -  2Q  •  K)  (*tt) 

vhere  a)  is  given  in  Eq*  (29). 


*Since,  in  general,  K  -  K!  is  not  parallel  to  K  or  K' . 
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In  order  to  obtain  an  order  of  magnitude  estimate  of  the  K  inte¬ 
grate  in  Eq.  (43)  we  take  Q  parallel  to  K.  With  this  we  find 


J  K2(  Q2  +  K2  -  2Q*K)  •*  3  ™ 


m 


By  substituting  Eq.  (44)  into  Eq.  (43)  and  perfoiming  the  Q  in¬ 
tegral  by  means  of*  the  delta  function  we  find 


r2  = 


2,31  Kmax  ( T  /m,  x2  5  f1  dx 

- w  (I/NM  u)  )  p  —* - 57T7o 

15( 2jr) 3  °  Ja  (x2  -  a2)1/2 

{x>-yjT7)\  ] 


(i5) 


where 

a2  =  2m  a)/p2  (46) 

The  lower  limit  on  the  integral  in  Eq.  (45)  canes  frcm  the  requirement 
that  K  be  real.  The  integration  over  x  is  elementary.  The  final  re¬ 
sult  is 


4m  \ax  (l/M  j.2  5(1  .  a2)l/2(i6/3  _  (l6a2/3)  +  a4)  (47) 

2  3  (2*)3  ° 


Putting  M  =  mass  Fe  atom  and  K  ~  100  R  we  find 

max 


r2  ~  3  x  io-11n 


or 


T2  =  X/r2  ~  10~  S6C 


m 


Next,  consider  the  process  shown  in  Fig.  3*  That  is,  double 
phonon  emission  by  two  successive  first-order  emissions;  each  of  which 
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is  characterized  by  an  amplitude  given  by  Eq.  (12)*  When  the  electron 
in  the  intermediate  state  is  nearly  free,  the  two  successive  anissions 
axe  essentially  independent  of  each  other  (i.e.,  this  case  is  essential¬ 
ly  multiple  single-phonon  emission,  and  should  not,  strictly  speaking, 
be  thought  of  as  double-phonon  emission).  In  this  case  ve  expect  the 
decay  rate  to  be  of  the  same  order  of  magnitude  as  r^. 

On  the  other  hand,  when  the  electron  in  the  intermediate  state 
is  highly  virtual  the  situation  is  quite  different.  For  example, 

Ein  "  Eint  ~  V  Where  Ein  iS  the  ener&  of  the  electron  ln 
initial  state,  and  is  the  energy  of  the  phonon  (which  is  in  a 

free  state)  and  electron  (which  is  in  a  virtual  state)  in  the  inter¬ 
mediate  state.  In  this  case  an  order  of  magnitude  estimate  for  this 
decay  rate  is  ~  r^/Ep  ~  10' 18  R.  This  rate  is  many  orders  of  magnitude 
less  than 

In  Section  IV  we  have  seen  that  for  steel  ~  10  ^  sec  ~  10). 

Thus  ve  conclude  that  double-phonon  emission  will  be  important  only 
when  ve  axe  dealing  with  a  very  pure  metal  (not  an  alloy,  since  for  a 
typical  alloy  the  unit  cell  contains  atoms  with  not  too  small  mass 
differences — compared,  say,  with  steel — and  hence  will  dominate). 

Lastly  ve  emphasize  that  the  calculation  °f  ^  Is/  very 

crude  since  ve  have  even  less  knowledge  of  the  second  derivative  of 
the  electron  potential  energy  than  we  have  of  its  first. 
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